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We propose to consider noniinear fluctuations in the theory of iiquid ^Fle deforming the commutation reiations 
between the generaiized coordinates and momenta. Generaiized coordinates are coefficients of density fluctu¬ 
ations of Bose particles. The deformation parameter takes into account the effects of three- and four-particie 
correiations in the behavior of a system. This parameter is defined from the experimental values of the ele¬ 
mentary excitation spectrum and the structure factor extrapoiated to T = 0 K. The numericai estimation of the 
ground state energy and the Bose condensate fraction is made. The eiementary excitation spectrum and the 
potential of interaction between the heiium atoms are recovered. 
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1. Introduction 

The method of collective variables is an effective approach to the study of an Af-particle Bose system. 
In this method, independent variables are the Fourier coefficients of the densi^ fluctuations of the 
particle (k is the wave vector). This method was originally proposed by Bohm It was suggested that 
from the whole set of an infinite number of variables pk, only DN variables are taken, where D is the 
space dimensionality. To achieve this, the domain of the wave-vector absolute values was restricted by 
some kc. 

An obvious imperfection of this approach is an ambiguous choice of variables. Bogoliubov and Zubarev 
@] proposed an approach where the set of values k is not limited. However, the transition from Carte¬ 
sian coordinates to collective variables pk can be provided by the weighting function. The hermitian 
Hamiltonian of the Bose liquid in the pk representation is written as a sum of Hamiltonians of an infinite 
number of non-interacting harmonic oscillators describing the oscillations of the Bose-liquid density plus 
a contribution from anharmonicities of these oscillations. Justifications of the feasibility of the method of 
collective variables are given in many articles using different Mproaches. Usually, the enharmonic con¬ 
tribution was considered using the perturbation operator IbHZdl. The respective results can be brought 
to the numerical calculations both for model systems and for strongly non-ideal systems like liquid ^He. 

A different approach for the study of the Bose system arose following the introduction of quan¬ 
tum spaces of minimal length. The simplest deformation is that of Kempf being quadratic in general¬ 
ized momentum The deformed commutation relations can be examined in multidimensional 

case. Such kind of deformation was applied to a wide variety of quantum mechanical problems, among 
which we distinguish the eigenvalue problem for D-dimensional isotropic harmonic oscdlator IZSll. th ree- 
dimensional Dirac oscdlator iz^l . (2-^l)-dimensional Dirac equation in a constant magnetic fiel(ij3 that 
were solved exactly. A composite system problem in the deformed space was also considered IZal^. 
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Some aspects of field theories were investigated in the deformed space, for example: the electromag¬ 
netic field a photoelectric phenomenon (ill], radiation and absorption of photons for deformed 

field jiil], the Casimir effect for the deformed field (i^l. An urgent problem is to find a possible gener¬ 
alization of arbitrary one-dimensional Heisenberg algebra with minimal length (or/and momentum) for 
multidimensional case jil]. Kempfs deformation can be applied not only to the operators of positions 
and momenta but also to collective variables. The idea to use a deformation of Poisson brackets to ex¬ 
plore a Bose system firstly appeared in the article fi^l . In paper lisll . the Bose particles were represented 
as a set of ^/-deformed harmonic oscillators. In paper (s^, there was used a representation of deformed 
creation and annihilation operators with a generalized four-parameter ^/-algebra. We note that the de¬ 
formed creation and annihilation operators associate with the non-linear /-oscillator operators. In paper 
fi^ . two-parameter deformed Bose gas model was proposed to find the correlation functions of the par¬ 
ticle. In a series of papers, the thermodynamics of ideal Bose and Fermi systems was studied. 

In paper (dll], the phenomenon of the Bose-Einstein condensation of the relativistic ideal Bose gas with 
deformed commutation relations for positions and momenta operators is investigated. 

In our paper, we associate the difficulties connected with enharmonic contribution with deformation 
of commutation relations between the positions and momenta. Thus, we replace the multimode Hamilto¬ 
nian of the Bose liquid in representation by a sum of single-mode Hamiltonians. We assume that this 
model successfully describes the properties of the Bose system when we make an appropriate choice of 
the deformation parameter. In this approach, we treat the collective variables as gerieralized coordinates. 
Thus, we attempt to take into account the enharmonic contribution described in with deformed Pois¬ 
son brackets quadratic in generalized momenta. In this paper, we use the deformation function which is 
quadratic in generalized coordinates. The correctness and efficiency of the proposed approach is shown 
for the analysis of Bose-systems with developed anharmonisms of the density fluctuations. Such a treat¬ 
ment of anharmonisms lowers the ground state energy if the deformation parameter is negative. It is 
interesting to note that the energy levels are bounded for such values of the deformation parameter. We 
make numerical estimations of the analytical results. 


2. Hamiltonian of Bose-liquid in the deformed space of collective vari¬ 
ables 


Let us consider a system of N spinless Bose particles of mass m and by the coordinates ri,...,rjv 
which move in the D-dimension space of the volume V. The Hamiltonian of the system reads 


iV A .2 




( 2 . 1 ) 


where the first term is the kinetic energy operator, pj - -\hVj is the momentum operator of the jth 
particle; the second term is the potential energy consisting of a sum of the particle interaction potentials 
®(|r; -r/). We introduce the collective coordinates representation which takes the form: 


N 


Pk- 


'/Nm 


k^iO. 


( 2 . 2 ) 


The Hamiltonian of the Bose liquid <2.1> can be rewritten in the following form US: 






^0 2m ( dpkdp_k 4 


1 1) N{N-1) N ^ 

+ 7PkP-k - 7 I + 777 — Vo + — Y Vk(PkP-k - 1) + Aff, (2.3) 


2V 


2V 


k^tO 


where 

Vk = J e“''^<t'(i?)dR 


(2.4) 
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is the Fourier image of the interaction potential. The operator AiT contains all anharmonic terms and in 
addition contains the term from the specific linear anharmonic by and quadratic by d/dp^: 

_ _ h^ikk') ^2 

= ZL -^Pk+k' 

k^ok'?to 2m VN 
(-)" 


dpkdpk' 


+ E 


H 


ki+...+k/2=0 


E ■■■ E :t— (fc? + --- + fc«)Pki...Pk„- 


(2.5) 


The introduced in (2.2> variables p^ are complex and can be represented as a linear combination of 
real variables: 


Pk 

Pk 


Pk-^Pk’ 

1 ^ 

Y coskr.', 

YNj^i ' 


1 " 


Since they are the complex conjugate value p* = p_k, that Is p^ = p^j^, p® = -p-k> and the variables with a 
particular index value k equal to (-k). Thus, Independent variables are pk when the wave vector k takes 
up values to the half-space of the wave vectors. Taking into account aU the remarks we have made, let us 
consider the harmonic part of the Hamiltonian of the Bose-liquid 12.3) as an infinite set of non-interacting 
harmonic oscillators with oscillation frequency cj^: 


E E' 

M=r,sk^tO 


IP. 


k,M 


2mk 


NIN- 1) ^ ( 

+ ^^vo-VI 
k^tO 


^ N 1 

1 4m ’ 


( 2 . 6 ) 


where the generalized momentum operator is conjugate to the generalized coordinate Qk,/j- In the 
coordinate pk-representations 

Qk,c=p^ Qk,s = p^ (2.7) 

Exphcit form of the generalized momentum operator is as follows: 

d 

PkM = -ifl^ -■ ( 2 . 8 ) 

dpk,fi 

The nteans that the sum takes the values k just from the haft-space of the domain. Comparing 
Hamiltonians (2.6) and <2.3) we have 

2m hk^ I ziv /h'^k^ 

We note that the operator Pk,^ has the dimension of action h, and Qk,^ is dimensionless. 

The Hamiltonian of the system <2.3) contains harmonic terms. We suggest that the anharmonicity AH 
can be taken into account by deformation of the commutation relations between generalized coordinates 
and momenta: 

Qk.^^Pk.^^ - Pk,^^Qk.^ = + PkQlf,), ( 2 . 10 ) 

where the dimensionless deformation parameter depends on the absolute value of the wave-vector. 
The operators and Pk'.fi' with different index commute. The deformation parameter for ^He can take 
negative values. We do not expect a full description of the properties of the Bose liquid Hamiltonian (2.6) 
with the condition (2.10) but assume that this model describes the behavior of the many-boson systems. 


3. Energy levels and wave functions of a Bose liquid 

We find solutions of the stationary Schrodinger equation from harmonic oscillators in the deformed 
space with Hamiltonian <2.6) . We suppose that the deformation is positive. For this purpose, we use the 
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canonically conjugated operators with the standard Heisenberg algebra: 


Pv .,^1 — Pk,fi ! 


Qk,^ — 


tan (4k,^ VW) 

7K 


(3.1) 


To solve the eigenvalue problem we use a representation of the deformed operators Hk,/j and Qk,/j (2.10) 
which express them in terms of canonically conjugate operators 4k,Pk.fi- The representation we have 
chosen reads: 

^k,/jPk,/j ~ Pk,fi^k,ii — ih. (3.2) 

Taking into account the representation (3.1) we rewrite Hamiltonian (2.6) in the form: 


E E' 


Pk,/i tan^ (4k,^ y/fe) 

2mk "^2 Pk 


N{N- 1 ) 

+ -Vo ■ 

2H 


^ [ 4m 


k^tO 


(3.3) 


The eigenvalues and the wave functions of the harmonic oscillator with deformed Heisenberg algebra 
were calculated in @]. Taking into consideration our notations (2.9) we find: 


E 


•••> ^k,c» ■ 


^^k,s> ••• 


E E 


2 m 


(«k,M + j] 


P=C,5k7t0 


1 + 


El]' 

2akj 


N(N-l) 

H-Vo 

2 V 


-I( 

k^tO 


2 1.2 




N 


I - 

y 4 m 2 V 


)• 


(3.4) 


here, quantum numbers «k,f; = 0,1,2,..., q - c,s. In the case of positive deformation parameter, the 
energy spectrum is infinite. The wave functions in the coordinate representation 4k,^ = Ik.fi, Pk,^i - 
-\Hdldqy_ ^ can be written as follows: 


V'...,file,^ 7k,fj. ■■■)- n V'nk./i Wk,^i) ■ 

k^o P=c,s 


(3.5) 


The ground-state wave function (Uk,^ = 0) reads: 


Vo(.qk,yt) 




r(v-ti) 


r(i/2)r(v-n/2) 


COS q, 


for 5:1, we have 


Vn{qk,yt) = P 


1/4 


r(v-i- n+ l)r(n-i-2v) 
n!r(l/2)r(v-t ?2-H/2)r(2v-t2n) 
d 


X f-h vtanq] f-h {v+ n - Ijtanq] cos''^” q, 

I dq ! \ da > 


dq 


here. 


= i + 2l,/, + (iL) 

2 Bi. \ (2ai.J 


2 Pk V ^ 2 q:j;( 

q — qk./i \lPkt ^ ~ ^k,ii ■ 

The wave functions are orthonormalized: 


71/(2 v/^) 

J" y^nPqk,ij)y^niqk,fP ^qk,fi — ^n',n ■ 

-nl{2s/Tk) 


(3.6) 


(3.7) 


(3.8) 
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The explicit form of the wave function of the first excited state is as follows: 

- cos'' q sin < 7 . 




2T{v + 2) 


r(i/2)r(v+i/2) 


The energy spectrum is quadratic in quantum numbers similarly to the theory of anharmonic oscil¬ 
lator in the case it takes into account the terms proportional to ~ lal^. 

Now we consider the case < 0. Then 


and we impose the requirement 




(3.9) 


(3.10) 


(3.11) 


The uncertainty relation obtained for the deformed algebra <3.9t leads to the fact that minimal uncer¬ 
tainty for the momentum operator is equal to zero. The new canonically conjugated variables (3.11 and a 
structure of Hamiltonian (3.31 with the changes -^\Pk\, tanC^k^^ '/Wic) tanh(4k,^i 'J\Pk\) whl be the 
same. Using the factorization method @| we can write the wave functions of the transformed Hamilto¬ 
nian for modes: 

/r(v+i/ 2 ) 1 


here, the brackets ((.. .)> denote an average: 

([■■■)) — J" y^oiqk.ii) (■ ■ ■) (^k,^)dt7k,fi • 


Vo{qk,ti) = \Pk\ 


r(l/2)r(v) cosh''^7’ 


when n>l-. 


Vn{qk,^i) 


\Pk\ 


1/4 


r(v-?i-f i/2)r(2v-2n-f 1) 
«!r(l/ 2 )r(v- n]r{2v- n+l) 


dq 


X-h V tanh -- V {v - n + \) tanh ^7 cosh 


dq 


q: 


here, 




2^ \Pk\ 


qk,ii 


\/m, 


\2ak> 
n = «k,/j ■ 


(3.12) 


(3.13) 


In case of deformed algebra (3.9) . the quantum numbers are limited by the number n<v and the energy 
levels are bounded from above. When v ^ n, the energy spectrum is continuous. The wave function is 
also normalized. We note that the integration variable ^7k,^i runs from 00 up to -00 


00 

/ 


Vniqk,fPVn(qk,^) ^qk,fi — 1 - 


(3.14) 


For example, we get for n-l: 


ViiqKfi) = \Pk\ 


1/4 


2 r(v-1-1/2) sinht/ 
r(l/2)r(v- 1) cosh''^ 


In the non-deformed case iPk = 0), we obtain the energy levels of the Bose liquid in the main approxima¬ 
tion @]. When in equations (3.61 and (3.121 . for the wave functions, the quantity pk—^Q and v ^ 00 , we 
obtain: 


cos'' q - 11 - 




\aklh 
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l/ifcl- 


cosh '' q - (1 + 
' -o'. 




\-ak/\Pk\ 


\lik\-0 




Taking into account the asymptotic formula for the gamma function r(v + a) ~ VSe when 

V ^ oo, we obtain the wave functions of the harmonic oscillator: 




= 


1/4 


1 


71 > 




^k,/i 


-7)^12 


d77 


(3.15) 


here, 77 =\/^, ak=mkWk/h. 


4. Ground-state energy 


If the quantum numbers in equation <3.4> are equal to zero, Uk,/i = 0, we obtain the ground-state 
energy Eq which can be written after transformations in the form; 


£0 = 


N{N- 


2V 


1 ) ^ fl^k^ 

—+E 

^ 8m ^ 


k^O 


k^to 


4m 


ak 



(4.1) 


The first two terms recover the ground-state energy in the Bogoliubov approximation Ql- The last term 
in equation <4.11 takes into account the anharmonicity and always leads to the negative values for the 
parameter Pk- The lowering of the ground-state energy of the liquid ^He can by achieved by a direct 
consideration of the enharmonic operator AH in the Hamiltonian (2.3) from perturbation theory EilE^ . 
The Fourier coefficient vo in the limit A; ^ 0 can be expressed via the speed of the first sound for liquid 
^He (when T = 0 K, the speed is c = 238.2 m/s). The definition of the speed of the first sound when T = 0 K, 


N 


d^Eo 

dN^ 


me 


allows us to obtain the equation: 

2 N 1 ^ h^k'^ (.al-lf 1 

mc^ = —Vo-E- 5 - „ . (4.2) 

4N^^4m a\ (1 + (/ifc/2afc)2)3/2 

Having taken into account the obtained equation (4.2) we rewrite the ground-state energy in the thermo¬ 
dynamic limit: 


£0 = 


Nmc^ 1 ^ h^k^ 


k9^:0 

1 ^ /z 2 fc 2 


-{ak-ir+- 


1 


+ — E 
16 ^ 


k9^:0 


2m ak 


-h--) 

ak’ 


1 -t 


_ h^k^ 

E t;—“ fc 

1^0 2 ??^ 

2,-3/2 

\2ak) 




1-tl —I + — 
^2ak’ 2ak 


(4.3) 


The found expressions can be used for the models with exact or perturbative solutions. To verify the va¬ 
lidity of the expression (4.3) . we might compare it with the results given by an exactly solvable model or a 
model which allows perturbative consideration. Such a comparison brings about an interesting question 
concerning the choice of deformation parameter pk- 


5. Elementary excitation spectrum 

Let us find an expression for an elementary excitation given by the wave vector q. Suppose that 
only the quantum number nq c = 1 and the other quantum numbers = 0 when k q, ^ c. From 
equation (3.4) we obtain: 

£...,0,nq,c = l,0,...;...,0,... = £...,0,...;...,0,nq,s=l,0,...;...,0,... -Eq + E^, 


33002-6 
























Theory of many-boson system with deformed algebra 


where the elementary excitation spectrum 


En — ^ n 

^ 2m ^ 




yiaq! 


If /3fc = 0, the latter equation leads to the Bogoliubov’s excitation spectrum 


(S.l) 


E 


B 

q 


2m 


i-q- 


(5.2) 


Equation (5.1) is an exact solution of the Schrddinger equation with Hamiltonian (3.3) . Here, we do not 
obtain the fading phenomenon of the elementary excitation because in Hamiltonian (3.3) there are no 
terms that describe the dissipation and the decay of the elementary excitation. The decay of the elemen¬ 
tary excitations in the liquid ^He can be explained by the fact that the spectrum of Bose-liquid ends at 
A;^3.6 A"i (Hi. 


6. Structure factor 

The structure factor of the system can be defined as an average of the square of the density fluctua¬ 
tions: 

Sfc = (|pkl">. (6.1) 

Now, we calculate the average (6.11 taking into consideration the fact that the collective variables are 
generalized ones and using the ground state wave functions (3.61 . Since the collective variables are gen¬ 
eralized coordinates from equation (2.71 . we find the average (6.11 on the wave functions of the 
ground-state (3.61 : 


Sk = 


7t/(2v^) 

Z<0k,^> = 7r^ f Vo(^?k,|i)tan2(y^pk,/ 

i=c.s Pk u=c,s J ^ 


^=c,s 


( 6 . 2 ) 


-7r/(2v^) 


r(v-f 1 ) 


PkTiv+l/2)Til/2) 


r(v-f 1 ) 


jr/2 

I 

-7tl2 

Jil2 


cos^^ qtan^ qdq 


PkTiv+ii2)rai2) 


J sin^pcos^^^ qdq. 


1 


-71/2 

This integral is reduced to the beta function: 

1 r(v-l/2) 

^'=“^r(v+i/2) “ Pk(y-i/2}' 

The parameter v should be taken from the relation lsTTl and after some transformations we obtain: 

1 




ak \/l + {Pkl2ak)^ 

When Pk - 0, we recover the well-know result by Bogoliubov and Zubarev ji|], 

Sfc= l/CTfc. 

The elementary excitation spectrum (5.11 reads: 

En^ — — + —^Pq. 

^ 2mSq 2m ^ 


(6.3) 


(6.4) 


(6.5) 
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The structure factor can be written as follows: 


2V SEo 

Sfc-1=— 

N 5vk 


here, F is the free energy of the system, which is equal to the ground-state energy £0 at T = 0 K. Having 
performed elementary calculations with equation ( 14.H we obtain the result <6.31 . Note that the struc¬ 
ture factor Sfc is an analytic function of the deformation parameter This fact can be seen when we 
calculated Sjc and make a parameter -\Pic\m the wave functions I I3.12> : 


Sk 


(X) 


00 „ 00 

2 r(v-tl/2) r tanh^t? _ 4 r(v-H/2) r sinh^^ 

I^Sfcl r(i/ 2 )r(v) J cosh2''t7‘^^“ I; 6 fci r(i/ 2 )r(v) j cosh^^^^ 

-<X3 ' 0 


( 6 . 6 ) 


This integral is reduced to the beta function: 

^ r(v-n/2) 1 

^ iPkir(v + 3/2) iPkl(v+l/2)' 

here, the parameter v is taken from the equation ( 13.13) . Having performed some simplifications we obtain 
the structure factor from the equation ( 16.3) . We note that condition )3.10) leading to the restrictions on 
the deformation parameter can be represented in the following form: 

IPkISk 

^ < 1. (6.7) 

2 


7. Potential and kinetic energy 


Having used the structure factor we rewrite the ground-state energy of the many-body system ( 14.1) : 


Nmc^ 1 y p 1 ■ 5 -^ 

£0 =-+ —E 

9 16 ^ 


k,tO ■ 




' k^tO 


1- —I +A£o, 

Ok) 


(7.1) 


here. 


1 h^k^ n iPkSk\MPkSk\ l^h^k^ 
A£o = —E- 2(Sl-V,+ \^-^\ -p-y- 

16^o2mSfc 2 2 4^^ 2m 




The mean value for the potential energy can be calculated for the ground state using equations dO 
After some transformations, we arrive at the relation for the mean value of the potential energy: 


(7.2) 

dO}. 


N{N-l) N _ p 

<®> = —TfT—vo-P —E Vfc((IPkl >-l) 

/V ^ ^ k^to 

N[N-Vi ^H^k^ p 


2V 


k^to 


nfc \/l-P(;Sfc/2aj;)2 


(7.3) 


Let us rewrite the average potential energy in terms of Sfc and speed of the first sound: 


($> = 


Nmc^ 1 ^ h^k^ 


16 


E 

k^O 


2mSk 


1)^-7 E 


2^2 


h'^k 


4kTo 


■(Sfc-1) 


1-^| + A(®>, 


(7.4) 


here. 


1 h^k^ p p 1 h^k^ IpkSk 

A(<1)) = — E- piiSk - lf{Sk -P 2) -p — E- 

32^0 2m 16k^„2mSfcl 2 


’ kptO ‘ 


(7.5) 
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The mean value of the kinetic energy is the average value of the operator I l2.6t in case = 0: 


E E 


/'=c,sk9to 


,h^k^ 

Am 


1 




(7.6) 


We calculate the first term in the square brackets on the ground state wave function from the equa¬ 
tion (3.6): 




r(v-ti) 


r(i/2)r(v-ti/2) 


nl{2y/h) 2 

\/^ / COS 


(7.7) 


= 


r(v-ti) 




T[\I2)T[V+\I2) ' 

v^r(v-i- 1 ) 




J 

cos i7k,u 




(7.8) 


r(i/2)r(v-n/2) 

( 

After simple transformations we get: 


7tl2 


sin'^ <7 cos^'' ^ qdq. 




L , PkSk? 


v-l /2 2Sk[ 2 I 


(7.9) 


here, the expression for v is taken from equation (3.7) and the structure factor from equation (6.3) . 
The second term in equation (7.6) is the structure factor by definition. Finally, we can write: 


<^^>=E 


I ak[\/T+{fij2a^+ Pkl2ak]^ 


k^O 


8 m 


1 


^/l + (Pk/2ak)'^ 


aic y/\ + {pij2a^ 


:-2 


Taking into consideration equation (6.3) we obtain: 


1 ^ (1 - Skf' 

(io^-E ^ 


4k^o 2 ??^ 


Sk 


MK). 


1 ^ h^k^ I PkSk] 
^ 2m ^ '' I 




4 ) 


(7.10) 


(7.11) 


(7.12) 


This result for the kinetic energy can be obtained as a derivative from the ground state energy with 
respect to the mass: 

dEo 

{K)^-m^. (7.13) 

dm 

The sum of expressions (7.4) and (7.11) is the total energy (7.1) . 

Now we find the potential of the interaction between the Bose particles: 


cl>(r) = 


(i^/ 


e''"vfc dk, 


(7.14) 


here, Vjt is the Fourier coefficient of the interaction potential between Bose particles represented as a 
function of the structure factor (6.3): 


Vk-- 


V h^k^ 
N Am 


E 1 

S? 4 


\^k ^ / 

Thus, the interaction potential between the helium atoms reads 


®(C) = 


1 




2ji^pr 2m 


CX. 

/ 


A:'^ sin kr 


1 


52 ^ 

\^k y 


dfc-i-A<I>(r), 


(7.1S) 


(7.16) 
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here, the contribution is A^ir) caused by the deformation of the commutation relations: 


A<h(r) = - 


1 

8ji^pr 2m 


OO 



sinfcr dfc. 


(7.17) 


8. Momentum distribution 


To find the average number of atoms with the momentum Mu, k 0 it is sufficient to calculate 
the variational derivative from the free energy of the system with respect to the free-particle spectrum 
12m. We note that the free energy of the system coincides with the ground state energy when T = 0 K 
[equation 14.3) by = 0], After simple calculations we get: 



Q:fc[\/l + (;6fc/2Q:fc)^ + j8fc/2Q:fc]^ ^ 

\/l + ()Sfc/2afc)2 


1 

Ufc \/l + {Pk/2ak)^ 



or after some transformation 



( 8 . 1 ) 


( 8 . 2 ) 


The expression for the kinetic energy can be represented as a function of average numbers of particles: 


E 

k^tO 


2m 


Nk- 


(8.3) 


The average number of atoms with the momentum hk can be obtained from the expression for the 
mean value of the kinetic energy <8.3) . To do so, we take into account the relation d7.10) . After simple 
calculations, we obtain the expression for the average numbers of atoms which coincides with the 
expression dS.l) . We estimate the relative number of atoms in the case of their momenta being equal to 
zero (Bose condensate): 


N Wk^o Sfc N 


(8.4) 


ANq 

N 


- 4^ E 

k^O 



PkSk ^ 


(8.5) 


At intermediate calculation, we assume that /3fc is of positive value. For numerical calculations of the 
deformation parameter it is negative. 


9. Deformation parameter. Numerical calculations 


For the numerical evaluation of the deformation parameter we proceed from the expression of the 
elementary excitation snectrui n d6.5) . The value of the structure factor and excitation spectrum are taken 
from experimental papers l4Si 14611 : 

h^ei2m~Tk' 


We have the values of the deformation parameter for a limited range of wave vectors because the elemen¬ 
tary excitation spectrum of the Bose liquid has the ultimate point of completion, and the experimental 
data of the structure factor are given up to 7.3 However, as fc ^ oo the elementary excitation spec¬ 
trum should be equal to the free-particle spectrum Ek h^k^l2m and the structure factor Sk 1. In 
figure [l] the deformation parameter Pk based on the experimental data for Sk and Ek is shown. We 
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Figure 1. Deformation parameter <9.H . 


shall model the deformation parameter by the function where a free parameter will be taken from the 
extrapolated data for the structure factor of the Bose liquid at T = 0 K 1^ : 

Pk = -Sk\Sk-lf- (9.2) 

The sign of the deformation parameter determines the behavior of the elementary excitation spectrum 
we have received (6.S). Therefore, Pk <0 because the experimental data for the elementary excitation 
spectrum of the liquid helium (triangles in figure [3] right panel) is lower than that of the theoretically 
calculated Feynman’s spectrum (circles in figurej^ right panel). 

The graph of the function (see ecruation l9.2) is shown in figure [2] (left panel). Note that this choice of 
the deformation parameter gives a correct behavior in the long-wavelength domain. Function 1 19.2t can 
be used to find the physical quantities in the limit of T ^ 0. The form of the curve for Pk (see figure [2) 
should be such that the elementary excitation spectrum (solid line in figure is reproduced for all the 
values of the wave vector. The point A: = 1 corresponds to the typical maximum on the curve of the 
elementary excitation spectrum (see figurerighf panel: triangles and circles). Thus, Pk (see figure[2) 
should have a clearly expressed minimum in this point. In figure[2](right panel) the dependence reflecting 
limitations imposed by the deformation parameter 16.7) is shown. 




Figure 2. (Left panel): Model of the deformation parameter )9.2> : (Right panel): The function = p^Skl2 
as a condition to limit the deformation parameter <9.21 . 

One can offer the model functions satisfying these considerations provided that contributions to the 
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basic physical quantities of the system reproduce the results in the post-RPA approximation in the theory 
of liquid helium-4. 

With the deformation parameter from <9.2> . numerical calculations of the obtained quantities per 
particle are made. In the thermodynamic limits (V ^ oo, N ^ oo, N/V - const), 



dk, 


the contributions of the physical quantities connected with the deformed commutation relations are es¬ 
timated. We do not expect a complete agreement between the properties of the Bose liquid within the 
proposed method and the perturbative results. However, we obtained the results that are consistent with 
those of the perturbation theory. The ground state energy with the deformation taken into consideration: 


N N~^ N ’ 


A£o 

N 


-1.89 K. 


(9.3) 


The numerical value of the first term corresponds to the ground state energy in the Bogoliubov approxi¬ 
mation, see (with p - 0.0219 and speed of the first sound c = 238.2 m/s): E^/N = -5.31 K. Thus, 

the fuU energy per particle in the deformation case is Eq/N - -7.2 K. The experimental result of the en¬ 
ergy is Eq/N = -7.13 K. Note that the linear corrections to the ground state energy over the deformation 
parameter give a leading contribution. This choice of the deformation parameter d9.2) offers an insignif¬ 
icant amendment to the value of the Bose condensate fraction obtained in the zeroth approximation: 


No _N^ ANo ANq 

"iv “ "aT ’ ~ir 


(9.4) 


In the zeroth approximation, N^/N gives a wrong result: N^IN - -0.31. 




Figure 3. (Left panel): The potential of the interaction between helium atoms. Circles correspond to the 
non-deformed case equation 17.161 with AO = 0. Solid line is the potential with models of the deforma¬ 
tion parameter equation <7.161 with AO 0; (Right panel): The elementary excitation spectrum. Circles 
denote Feynman’s spectrum; solid line is the spectrum of the deformed case; triangles correspond to the 
experimental data for the spectrum li^l . 

In our method, the anharmonic contributions from higher correlations in the Bose-liquid are taken 
into account within the single-mode approximation. The Bose condensate fraction taking into account 
the correction AA/o/A/in the deformed case [see equation <8.5) 1 is calculated in the approximation of the 
one sum over the wave vector. In this approximation, the value A/q/AT is negative. 

Nonlinear terms related to the deformation are quadratic in the expression for the interaction po¬ 
tential. This contribution vanishes, whereas the linear correction from the deformation parameter to the 
elementary excitation spectrum gives a significant contribution at q 2 A“^. This correction does not 
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contribute in the long-wavelength limit (figure!^. The behavior of the elementary excitation spectrum in 
the long-wavelength limit was solved in 14811 of the method of two-time temperature Green’s functions. At 
least a two-parametric deformation should be sought for to consider anharmonic terms AH in Hamilto¬ 
nian <2.31 in a proper way. This issue wUl be a subject of our future studies. 
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Teopifl 6araTo6o30HHoV ci/iCTeMi/i 3 AC(|)opMOBaHOio a/ire6poK) 
raM3eH6epra 

1.0. BaKapnyP, r.I. flaHO^KcP 

^ Kac|)eflpa reopemtHOi (fiianKi/i, /IbBiBCbKi/iCi HattioHaabHwii yHiBepcurer iMeni iBana (hpaHKa, Bya. 
flparoMaHOBa 12 , /IbBiB, 79005 , VKpaiHa 

^ /IbBiBCbKi/iCi HattiOHaabHuii yHiBepcMTer iMeni iBana (hpaHKa, npupoflHMti/iCi Koaeflx, Bya. TapHaBCbKoro 107 , 
/IbBiB, 79010 , VKpaiHa 


BanponoHoeaHO BpaxyBarn HeaiHiHHi c|)ayKTyaLtii b reopii piflKoro ^He, fletfiopMyKDHn KOMyrattiitHi cniBBiflHO- 
LueHHa Mi>K yaaraabHeHHMi/i KOopfli/iHaraMn ra iMnyabcaMi/i. B aKOcri yaaraabHeHi/ix KOopfli/iHar o 6 paHO Koec|)i- 
pien™ (fiayKTyapii rycTi/iHn 6 o 3 e-HacTHHOK. flapaMerp fletfiopMapii, mo spaxoBye Bnai/iB rpw- ra HOTnpnHacTHH- 
KOBi/ix Kopeatmiri Ha noBeflinxy 6 o 3 e-ci/iCTeM, o 6 paHO Bi/ixoflani/i 3 eKcnepi/iMeHtaabHi/ix 3 HaHeHb flat) cnexipa 
eaeMenrapHi/ix 36 yflxeHb ra eKcrpanoabOBaHi/ix eKcnepi/iMeHtaabHi/ix flani/ix crpyKrypHoro (fiaKTopa flo tom- 
neparypi/i T = 0 K. 3 MOfleabHi/iM napaMerpoM flec|)opMapii npoBefleno Hi/iceabHy opinxy eneprii ochobhoto 
crany ra KiabKOcri Boae-KOHfleHcary, BifliBopeno cnexip eaeMenrapHHX aOyflxenb ra norenpiaa B 3 aeMOflii Mix 
aroMaMH ^He. 

KaioMOBi cnoBa: iHecfiopMOBaHa anreSpa raHjaenSepra, cneKip e/ieMeHTapni/ix 36yfl>KeHb pi^Koro'^He, 
Boae-KOHfieHcaT 
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